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In this paper we build on the ideas presented in previous works for perfectly transferring a quantum 
state between opposite ends of a spin chain using a fixed Hamiltonian. While all previous studies 
have concentrated on nearest-neighbor couplings, we demonstrate how to incorporate additional 
terms in the Hamiltonian by solving an Inverse Eigenvalue Problem. We also explore issues relating 
to the choice of the eigenvalue spectrum of the Hamiltonian, such as the tolerance to errors and the 
rate of information transfer. 



I. INTRODUCTION 
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In a quantum computer, it will be necessary to perform gates between distant qubits, whereas the strength of 
interaction tends to reduce with distance, such that it is impractical to interact them directly. A typical response 
is simply to apply a series of SWAP gates to bring the qubits together so that they can interact. This, however, 
is cumbersome and risks introducing significant errors. Instead, it has been proposed that an ancillary device be 
introduced to act as a quantum wire. This wire would be a chain of qubits, with a fixed interaction, capable of 
transferring a quantum state from one end of the chain to the other. Following the initial investigation of wires Q , 
and subsequent demonstration that perfect quantum wires exist 0, Q , a large number of papers have been published 
about optimising the schemes over a variety ofparameters such as the robustness against errors, or a restricted ability 
to engineer the state (see, for example, 0,|EIE0)- Novel modifications of such chains have also been presented for 
the generation of entangled states or the application of unitary operations during the transfer || . The overhead of 
local SWAP gates is thus replaced by an engineering requirement. Such engineering can, however, be tested before 
the chain is used in a practical situation. 

None of these previous works have demonstrated perfect state transfer in a system that has realistic couplings, 
facilitated by dipole-dipole or Coulomb interactions, for example. Once such a coupling is introduced, the fidelity of 
such schemes is reduced below unity ,9j . Alternatively Q , two of these imperfect chains can be used in parallel to 
perfectly transfer the state. However, the trade-off is that there is no definite time at which the state has arrived. 
Instead, there is a measurement which will reveal whether the state has arrived or not. In this paper, we show how 
to adapt to arbitrary coupling schemes, hence pushing transfer schemes towards physical realisation. This is achieved 
by an iterative algorithm founded on the concept of Inverse Eigenvalue Problems (IEPs). The relation between 
IEPs and perfect state transfer has previously been noted in 0jE3- Given this relationship, it is also important to 
* -y understand why we might choose a particular set of eigenvalues, and how this affects the robustness of the chain when 
it experiences noise or other imperfections. We explore these issues in the second part of this paper. 
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?_j 1 II. INVERSE EIGENVALUE PROBLEM 

In order to make the connection to an IEP, we choose to make two assumptions. Firstly, by assuming the Hamilto- 
nian is spin preserving, [^2 <r Zl H] — 0, the problem is reduced to subspaces, and we can concentrate only on the first 
excitation subspace |3|]. By ensuring a single excitation is correctly transferred, a quantum state is also transferred 
because the state |00 . . . 0) is an eigenstate of the Hamiltonian. In the single excitation subspace, the basis states are 
denoted by |n), indicating the presence of the excitation on qubit n. 

Secondly, we shall assume the Hamiltonian is centrosymmctric (otherwise known as mirror symmetric). This means 
that for a chain of N qubits, the coupling between qubits i and j is the same as that between qubits N + 1 — i and 
N + 1 — j. Throughout this paper, we discuss chains of spins because these are the most efficient in terms of the 
number of qubits used. Further, following , we are assured that if the Hamiltonian can be mapped to a fermionic 
system (for example, the XY model can be mapped to a fermionic system by the Jordan- Wigner transformation), 
then by using a chain we get perfect state transfer in all excitation subspaces. 

The assumption about symmetry is useful because it ensures that the eigenvectors of the Hamiltonian are always 
symmetric or antisymmetric |l2j . It is well-known that a real tridiagonal matrix that exhibits these features has 
the additional property that when ordered with increasing eigenvalue, the eigenvectors are alternately symmetric or 
antisymmetric (see, for example, |l3f). There are no concrete statements of this form for the more generalised problem 
of non-nearest neighbor couplings. However, we can invoke the physical restraint that the coupling strengths should 
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drop off with distance, for example 1/r 3 . The symmetry of the Hamiltonian is maintained, so the eigenvectors are 
still symmetric and anti-symmetric. Only the eigenvalues are shifted. However, because the coupling strengths fall 
off with distance, they act as perturbations in comparison to the tridiagonal case. Provided these perturbations are 
smaller than the spacing between the original energies, then the Hamiltonian will still have the alternating structure 
of symmetry of the eigenvectors. This assures us that in this physically reasonable situation, a solution exists to the 
IEP, and that we know enough about the system to be able to generate perfect state transfer. 

Let us denote the eigenvectors of the Hamiltonian by |A n ), ordered such that the eigenvalue A„ > A n _i. The 
symmetry condition means that (i|A„) = (—1)" (N + 1 — i\X n ) ■ We are particularly interested in the case of i = 1. 
since this relates the initial state, |1), to the output state, \N). Starting with a single excitation at one end of the 
chain, we have |A„) = After a time t, the overlap of the evolved state and the target state, \N), is 

J2*ne- iKt (N\X n ) =J2<*ne- iKt (-l) n (1|A„) . 

n n 

Hence, if we apply the simple constraint on the eigenvalues that e _lA ™ to = (— 1)" for all n, we get perfect state transfer 
in time to. This has previously been observed for tridiagonal structures lldj . but we emphasise again that this 
applies to all centrosymmetric chains with coupling strengths that fall off with distance, or any other centrosymmetric 
network for which we have sufficient information about the ordering of the symmetry of the eigenvectors. 

The problem is now reduced to taking a desired eigenvalue spectrum, and a prescribed structure for the Hamiltonian, 
and solving for any free parameters that we might have (coupling strengths, site spacings, local magnetic fields etc.). 
Some classes of this problem are well-studied topics in the subject of IEPs j 1 21] - We now present a generalisation of 
the technique described in 0] , designed to cope with the arbitrary nature of H. 

Let us assume that we have a Hamiltonian H(a) which is represented by an N x N matrix in the first excitation 
subspace. This Hamiltonian depends on N parameters {a^}, ensuring that there are enough free parameters to be 
able to find a solution. Our desired eigenvalues are contained in the N x N diagonal matrix A [l^ . 

We start with a first estimate to a, a . The matrix H (a ) is diagonalised by Uq, 

H(J) = U A(l + eE )ul 

where E is a diagonal matrix which encapsulates the errors in the energies, and e is a small parameter. For our next 
guess, we will choose a vector a 1 — a + e 5a. Again, we can diagonalise the Hamiltonian, 

H(<?) = U 1 A(t + eE 1 )Ul. (f) 

We choose to parameterise XJ\ in terms of e, 

Ui = U {l+ieQ){l-ieQ)-\ 

where Q is a Hermitian matrix containing information on the change in eigenvectors. This parameterization ensures 
that Ui is unitary, and that UqU\ — > 1 as e — * 0. We can substitute this into Eqn. (JIJ, and expand in terms of e. The 
terms for e° cancel, so we choose to collect the terms for e . 



U Q = AE± - AE Q + 2i(QA - AQ) 



The aim of the iteration should be to choose 8a such that E± is minimised. Note that the diagonal elements of the 
final term, QA — AQ, are zero. Hence all the eigenvalue information is encapsulated by the diagonal elements of 
the equations, while changing eigenvectors only affects off-diagonal elements. Therefore, we select Ei to be the zero 
matrix, and rewrite the previous equation for just the diagonal elements. 

K.5a = e (2) 

e is a vector of the diagonal elements of — AEo, and the i th column of the matrix K is given by the diagonal elements 
of 
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FIG. 1: Example of 6 spins with local magnetic fields and coupled with a dipole-type coupling. 



The solution to Eqn. (01 is the vector 5a, which gives the correct eigenvalues to 0(e 2 ). Provided 5a is small in 

comparison to a , we can continue to iterate, squaring the error at each step. Hence, to achieve an accuracy of eo, 
we only need O(log(eo)) iterations. Since there are efficient algorithms for solving Eqn. (J2J, and because the matrices 
which we want to diagonalise are symmetric (hence there are efficient diagonalisation procedures, such as Householder 
reductions ^^)> the cost of each iteration scales polynomially (0(N 3 )) with the number of qubits in the chain [l5| . 
This means that we can solve for the required parameters with an efficient classical computation. 



A. Examples 

As an example, we could consider a system of 6 spins, arranged along a line, with separations and local magnetic 
fields as indicated in Fig. ^ As we have already seen, exact solution is possible if the couplings do not extend beyond 
nearest- neighbours. Instead, we shall now consider them to have a dipole interaction, so the potential ~ -^j. Hence 
the Hamiltonian in the first excitation subspace is of the form 
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Solving for the eigenvalues in this general form, and then equating these with the desired spectrum is certainly not 
an easy proposition. As we have seen, however, we can apply the algorithm described here, aiming for a spectrum of 
— |, —I, —\, \, §, §. The required values which we find are given in Tabled and are compared to those without the 
additional couplings. 



Parameter 


Exact Nearest- 


Dipole- Coupled 




Neighbor Solution 


Solution 


Bi 





0.491 


B 2 





-0.118 


B 3 





-0.373 


ri 


0.765 


0.967 


r-2 


0.707 


0.902 


rs 


0.693 


0.886 



TABLE I: Table of parameters for perfect state transfer for spin chain of Fig. Q 



One implicit assumption that we have made is that the Hamiltonian, H(a), is differentiable. There are some 
physical systems in which this might not be true. For example, we may be constrained to having to place spins 
on lattice sites of another material. Choosing which lattice sites to place the spins on is a discretized form of the 
problem, and is not covered in this formalism. The best that we can achieve is to allow some additional engineering, 
such as local magnetic fields, and tune these to give the closest match to a particular spectrum 0. We could envisage 
a variety of such systems in which we do not have control over a sufficient number of parameters. Instead of ./V 
simultaneous equations for N variables, we have N equations for m < N variables. We can solve for these variables in 
a least-squares sense, minimising the quantity AE±. Hence, while perfect state transfer might not be possible, we can 
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FIG. 2: Attempted perfect state transfer under various coupling scenarios. An overlap of 1 indicates perfect state transfer, 
(a) Uniform coupled chain, (b) Near-uniform couplings (within 1% of mean) with all spins coupled with a 1/r 3 interaction, 
without the ability to engineer local magnetic fields. High fidelity state transfer is possible, (c) Near-uniform couplings with 
ability to engineer positions of spins and local magnetic fields, enabling perfect state transfer, (d) The original perfect state 
transfer chain, which has a much wider peak and no periodic revivals at times t < to. 



maximise the fidelity of transfer. In Fig. [2] we have examined the case of N = 31, demonstrating that even without 
full control of a sufficient number of parameters, we can get higher fidelity than by simply assuming nearest-neighbor 
couplings. 

III. CHOICE OF SPECTRUM 

We now have an algorithm that takes a desired spectrum as input, and outputs the values of the parameters that 
we have access to in our Hamiltonian. It is therefore relevant to ask what spectrum we should choose. Unsurprisingly, 
the preferred spectrum is a trade-off between different properties. 

A. Coupling Variation 

One complaint that has been levelled at the construction of the original quantum wire is that the coupling 
strengths at the end of the wire are much smaller than those in the middle (by a factor of yN). This issue is 
relatively easy to correct by specifying a different spectrum, which closely matches that for the uniformly coupled 
spin chain. This selection of eigenvalues is still applicable to the ideas presented here. For example, we could take 
a chain of 31 spins, with uniform couplings and no magnetic field. We can solve for these eigenvalues, and truncate 
them to some precision (two decimal places in this example). We can then 'nudge' them slightly to give a spectrum 
suitable for perfect state transfer. The results we get require a variation in position of less than 1% about the mean, 
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and require local magnetic fields of the order of 10~ 2 . There are, however, trade-offs in terms of robustness against 
errors in the distances, and in terms of the time at which the arriving state is removed from the system. These 
problems are illustrated in Fig. Here we see that slight errors in the positions (i.e. the uniformly coupled chain) 
give far worse state transfer fidelity. We also see that the peak is very tightly confined with comparison to the original 
state transfer chain. 



B. Scaling of the State Transfer Time 

As presented in 0, the original state transfer chain has a fixed transfer time, to — n/2. However, there are two 
physical restrictions that we might consider applying to the chain which cause us to restate the transfer time. Firstly, 
we might choose to bound the maximum energy of the chain. In this case, we have to rescale the spectrum to fit 
within this bound, making the energy differences smaller. Since A max ~ N, when we rescale, we find that to ~ N. To 
minimise this scaling, we therefore want to choose a spectrum whose maximum eig envalue is O(N), and not 0(N 2 ), 
for example. Note that the uniformly coupled chain has a spectrum A& = 2 cos(fc7r / (N + 1)), which has a minimum 
spacing ~ for large N This minimum spacing corresponds to an upper bound to j- i.e. to ~ N 2 . 

The second approach is the most physically reasonable - we bound the maximum coupling strength in the system. 
For the original transfer chain, the maximum coupling strength is at the middle of the chain, ~ N. Hence, if we 
fix this to be constant by rescaling all of the coupling strengths, to ~ N. However, if we consider the chain with 
almost-uniform couplings, clearly there is no such variation. This is a very special case, and indicates that there is no 
direct correlation between general properties of a spectrum (such as E ~ ./V 2 ) and the coupling strengths. Since we 
can only choose the spectrum, it is very hard to apply this criterion. 

C. Timing Errors 

A further concern is the tolerance of the chain to a variety of errors. Perhaps the simplest case to consider is a 
timing error. If the perfect state transfer time is to, what is the fidelity of transfer, / = | (N\ e~ lHt |1) | 2 , at a time 
to + St? Clearly, it is preferable to have a broad peak, thereby maximising the tolerance to such errors. The fidelity 
of transfer at such a time is calculated from 

(JV] e -*^(*o+«) |i) = ^| a „| 2 (-l)™ e - lA "( t »+' 5 *) 

n 

W l-lStJ2K\ 2 (\n- Al) 

n 

(up to an irrelevant global phase) by taking the mod-square, 

c,2 n-l 
f = 1 - -^Y. l a "l 2 ( A » " Al ) E K| 2 (A„ - A m ) + 0(St 4 ). 

n m 

Given that we have ordered our eigenvalues such that A„ < A n _|_i, we can provide a simple lower bound for this 
quantity, 

/>l-^k| 2 (A n -A 1 f. 

n 

This bound is easily optimised by choosing a spectrum with minimum spread, and it is also clear that this will do a 
good job (although not perfect) of optimising the fidelity. The spectrum that fulfils this condition is just that which 
is used in the original chain [2(, with eigenvalues 0, ±2, ±4 etc. We refer to this spectrum as the SMS (spectrum of 
minimal spread). 

D. Manufacturing Errors 

One other class of errors that we need to consider are static errors, introduced at the manufacturing stage. If it is 
only possible to manufacture coupling strengths to within a certain tolerance, can we select a spectrum that minimises 
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the negative impact on the fidelity of transfer? We will consider, independently, the effects when this error propagates 
to errors in the eigenvalues and eigenvectors of the Hamiltonian. 

Let us first justify that errors in eigenvalues and eigenvectors are equally destructive. If a single eigenvalue is affected 
by an amount 8 (which is equivalent to an error in a coupling strength of order 8, independent of the eigenvalues), the 
fidelity is only affected by 0(8 2 ) - the first order effect only changes the phase of the incoming state, not the arrival 
probability. If an eigenvector is altered by errors on the chain, then it can still be written as a sum of a symmetric 
vector and an antisymmetric vector, by using a discrete form of 2g(x) — (g(x) + g(—x)) + (g(x) — g(—x)). Of course, 
we cannot only change one eigenvector, since we need to maintain an orthogonal set of vectors. We can test what 
happens when 

|Ai) -> v / l-^ 2 |Ai)+5|A 2 ) (3) 

|A 2 ) -» *|Ai>- Vi-£ 2 |a 2 >, 

which serves to indicate the general properties of the errors. In particular, we can calculate that the fidelity of transfer 
is still only affected to 0(S 2 ). Hence, errors in the eigenvectors are no more destructive than errors in the eigenvalues. 

Recall that the condition we imposed upon the eigenvalues was e~ ™*° = (—1)", which means that if there is an 
error in A„, we are concerned with how close \ n to/Tr is to an integer. Hence if A„ is of the form A'f eal + Sir /to, then 
5, the inaccuracy, must be small in comparison to 1, not A™ . When we quote the accuracy to which an eigenvalue 
must be engineered, we consider it as a fraction of A max . Hence, this fraction is largest for the spectrum with smallest 
maximum eigenvalue i.e. the SMS. 

To consider the effect of errors on the eigenvectors, and how to counter them is a more difficult problem because 
we do not have any direct control - our eigenvectors are determined by the choice of spectrum and the Hamiltonian. 
However, we can make some progress. Let us introduce the N x N symmetry operator, S, defined by 

Si J = S^jf+l—j- 

This operator commutes with the centrosymmetric Hamiltonian (this is how we prove the symmetry of the eigenvectors 
of H), and hence, should be effective at detecting any symmetry that is broken by the effects of noise, by calculating 
[H, S]. Let us, again, consider the effect described in Eqn. ©. 

[H, S] = 2SVl-S 2 (X 1 - A 2 ) (|A 1 )(A 2 | - |A 2 )(A 1 |) 

This tells us that a fault in the coupling strength J' n = J n + 5 will translate into an error in the eigenvectors of 
approximately 5/ (Ai — A 2 ). Minimising this error is achieved by maximising the energy difference Ai — A 2 . This would 
appear to counter the effect that we just observed for the eigenvalues. However, we note that in the general case, 
where we have to consider all of the eigenvectors being affected, the most significant error is related to the smallest 
energy difference which, in turn, is determined by the state transfer time, which we are holding fixed. Hence, the 
choice of spectrum has limited effect on how the eigenvectors contribute to a reduction in the fidelity. It is more 
important to optimise for errors in the eigenvalues. 

E. Rate of Information Transfer 

So far, we have considered only one simple protocol for transferring states with our chain; we place states on the 
chain at some time, t = 0, and either remove them from the chain, or interact them with another qubit (leaving them 
on the chain) at t = t$. However, the rate of transfer suffers at this stage because if we bound the energy of the chain, 
the state transfer time increases with N (in the case of the SMS). If we restrict ourselves to only being able to access 
a fixed number of qubits at each end of the chain, then the transfer rate ~ -L. Is there something better that we can 
do, perhaps allowing some fixed reduction in fidelity? We restrict ourselves, without loss of generality, to being able 
to access only a single qubit at each end of the chain. 

Consider placing an excitation |1) onto a chain with a state transfer Hamiltonian, H, at a time t — 0. At some time 
td < to, we shall discard the information on qubit 1, and replace it with a new quantum state. We need to evaluate 
the fidelity of arrival of the initial state when the replacement we make is to place the first qubit in the |0) or |1) 
state. The first period of evolution gives us 

N 
n-l 
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where (3 n {t d ) is the amplitude of the state at time td on qubit n. If the first qubit then gets indiscriminately reset to 
the |0) state, the resulting density matrix is 

p(U) = \Pi(t d )\ 2 |0) (0| + (e-* m * |1) - (3i{t d ) |1))«1| e im « - R(t d ) (1|). 

After the further evolution for a time to — td, we expect to be able to remove the state from qubit N. The state at 
this time is 

p(t Q ) = \(3i(t d )\ 2 |0) <0| + (\N) - MtdK 1 "^-^ \1))((N\ - p* x {U) (1| e^-**)). 
Given that the chain is centrosymmetric, Pn^q — td) = (3i{td)- Hence, we find that the fidelity of transfer is 

io = (l-|/Mi d )| 2 ) 2 . 

Similar arguments lead to a calculation of the fidelity of transfer when we reset to the |1) state. This is complicated 
compared to the previous argument because we now enter the second excitation subspace. The structure of higher 
excitation subspaces can be described by the eigenvectors through the Slater Determinant (see, for example, ^3). An 
alternative approach is to make use of the wedge product, whose application to these systems was described in |lfij . 
With a probability (1 - \(3i{t d )\ 2 ), we introduce a second excitation into the system, and the state is described by 

|i)A( e -^|i)-/M^)|i». 

The properties of the wedge product automatically mean that it is impossible to get two excitations on a single 
spin (|1) A |1) = 0), and the evolution of the state is described by applying the single excitation evolution to each 
component, 

-» ( e - iH( *>- td) A \N) 

N-l 

= J2^(t -t d )\n,N). (4) 

71=1 

We conclude that if the reset was made, the state definitely arrives. The other possibility, of having a |1), occurred 
with probability |/?i(id)| 2 and after a time to — td gives an \N) with probability |/3i(t d )| 2 . Hence 

fi = i-mtd)\ 2 + \Mt d )\ 4 
= h + \pi{t d )\ 2 

With this formalism, we can clearly calculate similar transfer fidelities in all excitation subspaces, and we find that 
the minimum fidelity is /o, provided /3i(i) is a decreasing function for t < to- By symmetry, the fidelity of arrival the 
newly added state at a time td + to on qubit N is the same, if the first state is removed at time to- These fidelities only 
depend on |/3i(id)| 2 , which is the fidelity with which the initial state is still on the first qubit at time td- Therefore, 
to maximise the rate of information transfer, we select the chain with the thinnest peak, the opposite of Sec. IIII CI 
but this already assumes that we have control of the timings to within the width of the peak. Fig. [2] also serves to 
illustrate another point that we have to be careful of. For the original state transfer chain, f3\{td) is a decreasing 
function for times td < to, and hence the effect if we also add an extra state at 2td, as well as at td, is small, to 
the point that it becomes insignificant. However, other spectra can give periodic revivals on the first spin during the 
transfer. We have to ensure that we avoid all of these, which potentially reduces the rate of information transfer. 

For the SMS, this protocol enhances the transfer rate to ~ -7=, for some fixed reduction in fidelity, e. This is an 
enhancement over Jl7| , where it was demonstrated how a uniformly coupled ring, in combination with Gaussian wave 
packets (to minimise dispersion), could be used to give a transfer rate, which was lower bounded by A. One point to 
note is that the process of removing the state at t = to (Eqn. (0}) means that only a single excitation is left on the 
chain. Hence, there are only ever 0(y/~N) excitations on a chain, and we never saturate the number of excitations on 
the chain. 



IV. CONCLUSIONS 

In conclusion, we have demonstrated that perfect state transfer is possible in the presence of next-nearest-neighbor 
couplings by presenting an algorithm that correctly calculates the couplings for any specified system. If sufficient 
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free parameters are not available, the formalism presented here is easily adapted to find the optimal solution in a 
least-squares sense. This technique may generally be useful when, as a first approximation, people have considered 
the restriction to nearest-neighbor couplings |a], and would subsequently like to extend their networks to include 
more realistic Hamiltonians, such as those met in physical implementations. We have also discussed some of the 
issues relating to what spectrum should be chosen for the state transfer, demonstrating that the spectrum originally 
proposed in [2| is close to optimal is terms of robustness against a range of errors. 

The author is supported by EPSRC, and the Denman Baynes Studentship at Clare College, Cambridge. He would 
like to thank Sonia Schirmer for useful discussions. 
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